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Abstract. Methods and properties regarding the linear perturbations are dis- 
cussed for some spatially closed (vacuum) solutions of Einstein's equation. The 
main focus is on two kinds of spatially locally homogeneous solution; one is 
the Bianchi III (Thurston's H 2 X R.) type, while the other is the Bianchi II 
(Thurston's Nil) type. With a brief summary of previous results on the Bianchi 
III perturbations, asymptotic solutions for the gauge-invariant variables for the 
Bianchi III are shown, with which (in)stability of the background solution is 
also examined. The issue of linear stability for a Bianchi II solution is still an 
open problem. To approach it, appropriate eigenfunctions are presented for 
an explicitly compactificd Bianchi II manifold and based on that, some field 
equations on the Bianchi II background spacetime are studied. Differences be- 
tween perturbation analyses for Bianchi class B (to which Bianchi III belongs) 
and class A (to which Bianchi II belongs) are stressed for an intention to be 
helpful for applications to other models. 



1. Introduction 

A classification of three dimensional homogeneous (Riemannian) manifolds is 
known in relativity as the Bianchi classification, which is a classification of the 
three dimensional simply transitive Lie groups by isomorphism. A four dimensional 
spacetime {X x R, g a b) is said to be spatially homogeneous if a Bianchi group G acts 
spacelike and makes the metric g a b invariant. Cosmological models based on such a 
spacetime structure, especially with isotropy, have long been a basis for cosmology 
and (in case of being anisotropic) have recently been subjects for more mathematical 
study in relativity (e.g. |WE| ). Departure from a homogeneity however belongs to a 
very recent development in the field (except for linear perturbations of rather special 
homogeneous solutions such as those with isotropy or ones that are not isotropic 
but can become so in a limit). In particular, whether the solution is stable or not 
is unknown for almost all such solutions. 
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This article discusses linear perturbations of two Bianchi type solutions; one is 
Bianchi III, the other is Bianchi II. The spatial manifold of each solution is assumed 
to be closed (meaning 'compact without boundary'), introducing an appropriate 
compactihcation. (As a result, the solution becomes locally homogeneous.) One 
of the common features of these two models is that both of the spatial manifolds 
belong to so-called Yamabe type —1, which is of special interest in relativity FM . 
One of the biggest differences is that while Bianchi III belongs to so-called class B, 
Bianchi II belongs to class A. (The class A or B is defined with respect to whether 
the trace of the structure constant of a Bianchi group vanishes or not. EM ) In 
fact, we will see some of the largest differences between Bianchi III and II amount 
to those between class A and B. Therefore, in particular, the basic methods for 
Bianchi II will be applicable to other class A Bianchi models, too. (Difficulties in 
Bianchi II will also be common for all Bianchi A models.) 

We deal with the Bianchi III perturbations first. In this case, we have al- 
most complete results. The technical details of fundamental analysis of Bianchi 
III vacuum perturbations, including the derivation of the wave equations for the 
gauge-invariant variables, have already been given in |TMY| and is also briefly 
summarized in subsections 12 . 1 1 to 12 . 31 of the present article. The asymptotic stabil- 
ity of the perturbations has been discussed in |YT| in a framework of dynamical 
system approach, where one does not have to appeal to an exact background solu- 
tion. In the present article we discuss the asymptotic stability from a different point 
of view in subsections 12 . 41 and 12.51 using the exact background solution. The issue 
of stability of the Bianchi II solution is on the other hand still an open question. In 
this article we discuss properties of some simpler field equations on that background 
in the subsequent section, stressing differences from the Bianchi III background. 

2. Bianchi III vacuum perturbations 

2.1. The background solution. First, let us describe the background ge- 
ometry briefly. 

The Bianchi III algebra is generated by three vectors (J = 1 ~ 3) such that 
the commutation relation is given by 

(2.i) [a, 6] = a, [6,6] = K 3 , a] = o. 

Using coordinates, they can be represented as 

„s d . d d , d 

2-2 fr = -fi-, ^ = !C «- + y«-' & = 

Ox Ox Oy Oz 

We will be able to think of these vectors as Killing vectors on our background 
spacetime. One can also consider invariant vectors Xl ancl 1-forms a 1 defined by 

(2.3) C^XJ = 0, and C^a J = 0, 

where denotes the Lie derivative with respect to the vector £/. Suppose xi ancl 
a 1 are dual to each other; {xii<^ J ) = S'J. Then, using the coordinates they are 
represented as 



(2.4) 
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For convenience we define a standard 2-dimensional hyperbolic metric 

(2.5) h = a 1 <g> cr 1 + a 2 <g> ct 2 , 
and a standard 1-dimensional Euclid metric 

(2.6) l = a 3 ®a 3 . 

Then, the (universal covering) vacuum solution we consider is given by 

(2.7) g = -N(t) 2 dt 2 + qi (t)h + q 2 (t% 

where the two scale functions qi(t) and q 2 (t) and the lapse function N(t) 2 are 
explicitly given by 

(2.8) qi (t) = (t + k) 2 , q 2 (t) = (t-k)/(t + k), N(t) 2 = l/q 2 (t), 

where k is a real parameter. 

This type of solution is called locally rotationally symmetric (LRS), because, 
due to the hyperbolic metric h the spatial metric q = q±h + q 2 l has a fourth Killing 
vector (with degenerate points). We consider only this type of metric, because 
otherwise one cannot compactify the spatial manifold. 

We can express the compactified solution as follows: 

(2.9) (MxR,g) = (R 4 ,g)/r, 

where the spatial manifold M is supposed to be the direct product S 9 x S 1 , where 
S g (g > 1) is a higher genus surface with genus g while is S 1 the circle. For 
convenience, we call E s the base and S 1 the fiber. T C Isom(R 4 ,g) is a discrete 
subgroup of the isometry group of (R 4 , g), which acts on (R 4 , g) from the left. We 
assume that T is of orthogonal tvpe |TMY| . which means that (the action of) T 
preserves the foliation by H 2 of the universal covering manifold. In this case one 
can write the spacetime metric g for the spatially compactified manifold in almost 
the same form as g: 

(2.10) g = -N{tfdt 2 + q x {t)h + q 2 (t)l, 

where h is the hyperbolic metric on S 9 induced from the universal cover metric q, 
while I is the Euclid metric on S 1 induced from the same metric. 

2.2. Eigenfunctions, vectors, and symmetric tensors. We need to find 
complete sets of eigenfunctions, (eigen)vectors, and (eigen) symmetric tensors for 
the Laplacian so that we can mode-expand the field variables on our background. 
Thanks to the fact that our spatial manifold has the orthogonal direct product 
property, we can reduce this task to those for the base and fiber. Namely, once one 
finds eigenfunctions, vectors, and symmetric tensors on both the fiber and the base, 
the remaining task is just to make appropriate products. For details, see |TM Y| . 
Let us in the following recall some basic definitions. 

The most fundamental are the eigenfunctions on the fiber and the base. Let 
c m be a real eigenfunction on the fiber, corresponding to eigenvalue m; 

(2.11) (X3) 2 c m = -m 2 c m - 

We call m the fiber eigenvalue. (x3 — d z is the derivative operator along the fiber.) 
We also define another real eigenfunction c m by 

(2.12) X3c m = -mc m , \3C~m = mc m . 
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Similarly, let S\ be a real eigenfunction on the base, corresponding to eigenvalue 
A; 

(2.13) A h S x = -X 2 S X . 

We call A the base eigenvalue. (A^ is the Laplacian with respect to the metric h.) 
The mode function on the spatial manifold M — E g x S 1 is given by the product 

Sm , A = C m S\ , 

(2.14) A g 5 m , A = -(Qi 1 ^ 2 + q^m^S^x, 

where q — q\h + is the spatial part of the metric g. 

The mode vectors and symmetric tensors on (M, q) are categorized, according 
to those for the base (E ff ,/i), into four kinds; the even type, the odd type, the 
harmonic type, and the transverse-traceless (tt) type. The meaning of this is as 
follows. 

The even vectors (or 1-forms) S a on the base are equivalent to the exact (gradi- 
ent) 1-forms on the base, S a = D a S, where D a is the covariant derivative operator 
on the base. S is generally a function on the base, but it is supposed to be the 
mode function S = S\ when thinking of S a as a mode vector on the base. Similarly, 
the odd 1-forms V a are equivalent to the (Hodge-)dual exact 1-forms, V a — e a b D b S, 
where e ab — 2a 1 [ a a 2 b ^ is the area 2-form on the base. (To raise indices, use h ah , the 
inverse of h a b-) The harmonic vectors U a are the Hodge harmonic vectors defined 
by dtl = SU = 0, or equivalently D[ a U b ] = D a U a = 0. The number of the inde- 
pendent harmonic one-forms is equivalent to the lst-Betti number of the surface, 
which is given by 2g in case of the higher genus surface. The Hodge-decomposition 
tells us that an arbitrary 1-form on the base can be uniquely decomposed using 
these even, odd, and harmonic (mode) 1-forms. 

The even symmetric tensors are symmetric tensors made from the even vectors 
or the metric h a b and the mode functions. (The mode functions themselves are 
supposed to be of even type.) We define the even trace part by Sh a b, and the even 
traceless part S a b as the traceless part of the gradient of S a . Similarly, we define 
the odd symmetric tensors V a b and the harmonic symmetric tensors U a b by the 
(symmetric) gradient of V a and U a , respectively. (They are automatically traceless.) 
From York's decomposition jYOj , we know that it is necessary for completeness to 
consider the transverse-traceless (tt) tensors W a b, which are defined by D b W a b — 
h ab W a b = 0. The number of independent tt tensors on the higher genus hyperbolic 
surface is the same as the dimension of the so-called Teichmuller space, which is 
equal to 6g — 6. 

The symmetric mode tensors on the spatial manifold M have, as a result, nine 
kinds, which are denoted as (Ei) a b (i = 1, • • • ,9). Four of them are of the even 
kind, which are defined as {Ei) ab = c m Sh ab , (E 2 ) a b = c m S ab , (E 3 ) ab = c m Sl a b, 
and (E 4 ) a b — 2c m S( a a 3 b y (We omit the subscript A like S = S\.) Similarly, 
the odd symmetric mode tensors are defined as (E^ab = c m V a b, and (Eg) a b — 
2c m V( a a 3 b y The harmonic symmetric mode tensors are defined as {E-j) ab = c m U ab , 
and (E$) a b = 2c m U( a a 3 b y Finally, the tt symmetric mode tensors are defined as 
(Eg) ab = c m W a b- These are orthogonal to each other with respect to the standard 
L 2 -norm. (The mode vectors on M are also defined in a similar way.) 
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2.3. Vacuum perturbations and wave equations. Using the symmetric 
mode tensors denned above, we can expand the first variation Sq a b of the spatial 
metric q a t as follows. 

(2.15) ^ o6 = ^7 i ( J Bi)a6, 

where the sum is taken over i = 1, ••• ,9, as well as all possible eigenvalues A 
and m. The coefficients 7* = j l (t) are functions of time that are supposed to be 
perturbation variables. 

An important fundamental result that can be confirmed by direct computations 
is about how the perturbations decouple, which is stated as follows. For given base 
eigenvalue A and the fiber eigenvalue ra, the set of the even variables, the odd set, 
the harmonic set and tt set evolve independently from each other. Note that this 
statement says two kinds of decoupling; first of all, different eigenmodes decouple 
from each other, and the four kinds of set of variables decouple from each other. 
Due to this decoupling property, we can deal with the perturbations one by one for 
each one of the four kinds for each mode. 

The perturbation variables change their form by an infinitesimal diffcomor- 
phism acting on the manifold. Certain combinations of the variables however re- 
main invariant. Such a combination is known as the gauge invariant variable. In 
general, one of the central issues in perturbation analysis is to find these gauge 
invariant variables and determine the dynamics of them. Our choice of the gauge- 
invariant variables is as follows: 

-AX7 1 - A 27 2 + 7 3 - (2m/A) 7 4 , 

0VV2) 7 5 +7 6 , 
(m/2) 7 7 + 7 8 , 

7 9 , 

where Ai = q 2 /qi, A 2 = (A/^2(A 2 + 2))(Ai + 2m 2 /A 2 ), and v = m/VA 2 + 2. The 
subscripts E, O, H, and T are attached to express it is for the even, odd, harmonic 
or TT perturbation, respectively. 

Getting wave equations (equivalent to the linearized Einstein equation) is a 
substantial part of the work in the study of perturbations. While we show only the 
results of it in the following, it is worth mentioning that the actual computations 
were done in the Hamiltonian formalism, using so-called the method of generating 
function to find a desired canonical transformation that is necessary to split the 
system into gauge-dependent and independent parts. After lengthy computations 
TMY , we get the wave equations for the gauge invariant variables, which are 
second order ODEs of the form 

(2.17) Q A + /a(*)Qa + <7a(*)Qa = 0. (A = E, O, H, or T) 

Each coefficient function /a(£) or <?a(£) are given as follows: 
(i) For the even perturbations: 

f(t -2k) X\ , t\ , 1 Y 



(2.16) 



Qe 
Qo 
Qh 
Qt 
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where 

X = 8m 2 (t - 2k)t+ 2 + 2 A 2 m 2 (2t - k) t + 2 + 2\ 2 (2t- 3 k) + A 4 (2 1 - k)t/t+, 
Y = 16m 4 t+ 4 -4km 2 t+ [A 2 (t-4k)/t- -8] +2 A: 2 A 4 (2t + k)/(t-t 2 + ), 
Z = 4 to 4 t + 6 +8m 2 L t + 3 + 4 A 2 m 2 1 t+ 3 + 2 A 2 t_ i+ + A 4 i 2 . 

(ii) For the odd and harmonic perturbations (A ^ for the odd, A = for the 
harmonic) 

2^ 2 (t- 2fc)t + - u 2(t-3fc) 4ft - 2fc) 

(2 ' 19) /o = ' - 9o = (A +2) ^ + ^T 

where u = f 2 ^ + t_t^ 4 . 

(iii) For the tt perturbations 

dirt t 2(<-2fc) 2*+ | 2(t - 3fc) 

(2.20) /x = ____, ffT=TO _ + ___. 

2.4. Asymptotic solutions. Let us find future asymptotic solutions to dis- 
cuss the future stability. A lot is known about how we can find asymptotic solutions 
for linear ODEs. The following simple result, taken from standard texts, is sufficient 
for our purpose. 

Theorem 2.1 (e.g., |CH| ). Consider the following second order ODE 

(2.21) X" + f(s)X = 0, 
where f(s) is a function of s that approaches a constant 

(2.22) lim f(s) = C = constant. 

(Primes 1 stand for d/ds.) If C ^ and 

(2.23) / \f{s) - C\ ds < oo, 



then the equation $2.21\) has a set of fundamental solutions 
when C > 0, or / e ±V^r s (l -|- o(l))| when C < 0. If C = and 



(2.24) / s\f(s)\ds <oo, 

then the equation 1^2. 21}) has a set of fundamental solutions {1 + o(l), s(l + o(l))}. 
With the aid of this theorem we can prove the following. 

Theorem 2.2 (Asymptotic solutions for the generic case). Assume that the 
fiber eigenvalue m does not vanish, m^O. Then, the wave equations \2.1 r l\j for the 
gauge-invariant variables ]2.1b)) possess the following fundamental solutions 

Q A (t) = {te ±i l« a+akh *Q (l + o(l))}. 
(A = E,0,H, or T) 
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PROOF. Define a new time coordinate s by 

(2.26) *_,(„. 1+ » +< x£). 
Rewrite the wave equation l|2.17[l in terms of s. Then, by taking 

(2.27) aA =* e -iJfA*t 

and putting Qa = o-aXa, one has an alternative equation X A + Wa(s)X a = 0. 
It is easy to check that J°° \ W A (s) — m 2 | ds < oo. Therefore, from the previous 
theorem, Xa has the fundamental solutions e ±ms (l + o(l)), from which the claim 
follows. □ 

Similarly, we can prove the following for the m = case. This case is called the 
[/(l)-symmetric case, because this kind of perturbation keeps the £7(l)-symmetry 
of the background along the S 1 fibers. 

Theorem 2.3 (Asymptotic solutions for the ?7(l)-symmetric case). Assume 
m = 0. Define p\ = A 2 — 1/4, for convenience. Then, the wave equations {2.17}) for 
the gauge-invariant variables $2.16)) possess the following fundamental solutions: 

For the even perturbations: 

f j_I_ e ±VP^(l + (i))} (A 2 >I) 
(2-28) Q B (t) = J + o(l)), ^(1 + o(l))} (A 2 = i) 

[{t-Wi^l (1 + (1))}. (A 2 <±) 

For the odd perturbations: 

f{Vte ± VpTi°g* (i + (i))} (A 2 >i) 

(2.29) Qo(*)= \ {V*(l + o(l)),V^logt(l + o(l))} (A 2 = i) 

[{t^VM (l + o(l))}. (A 2 <±) 

Proof. Define a new time coordinate s by 

(2.30) | _„».(«,. J +0(1). 
Rewrite the wave equation l|2.17|l in terms of s. Then, by taking 

(2.31) a A = ^=e~^ fAd \ 

and putting Qa = oaAa, one has an alternative equation X A + Wa(s)Xa = 0. 
(A=E or O.) It is straightforward to check that 

f°° 1 

(2.32) / \W A (s)-px\ds< oo (A 2 f -) 

and 

f°° 1 

(2.33) / s\W A (s)\ds < oo. (A 2 = -) 

Theorem 12 . II therefore applies, and the claim follows after replacing the time coor- 
dinate back to t. □ 
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Remark 2.4. As seen from the last two proofs, the key point to find asymptotic 
solutions is the choice of a good new time variable (like the choice of a(t) or <j(°'(t) 
in the proofs.). 

Interestingly enough, in the generic case the asymptotic solutions are common 
for all kinds of perturbations. Moreover, the asymptotic solution does not depend 
on the base eigenvalue A. In this sense the asymptotic solutions are very universal. 

On the other hand, for the £/(l)-symmetric case, for each kind of perturbations 
the behavior of the asymptotic solutions is divided into three cases, depending upon 
the value of A. 

2.5. Stability issue. Although all the asymptotic solutions above show that 
our gauge-invariant perturbation variables are growing in time, one cannot say 
anything about stability of the solution from these results themselves. This is 
because they do not take into account the fact that the background solution is 
expanding. Because of the anisotropy of the background solution, it is a very 
subtle question how one can subtract this expansion effect. If the background was 
isotropic there exists a natural way to normalize the perturbation variables, because 
it has only one scale factor, but in an anisotropic case it is not clear how to find a 
right way of normalizing the variables, especially in a way such that the result is 
gauge invariant. 

Nevertheless, one possible way may be to use the zero mode solution to nor- 
malize the gauge invariant variables. By a zero mode solution we mean a solution 
with the vanishing eigenvalues. Such a solution represents a perturbation from a 
locally homogeneous solution to another locally homogeneous solution, which is an 
effect we are not interested in. So, it may make sense to use the zero mode solution 
to do a normalization. Moreover, it is apparently gauge-invariant. (This scheme is 
essentially equivalent to the one discussed in |YT| .) 

All the zero mode (vacuum) solutions are already given in TMY , from which 
one obtains the following growing rates for those solutions: 

(2.34) Q^=0(l), Q$P=0{t), Q { ° } =0(t 2 ). 

(Although there does not exist a zero mode for the odd perturbations, we identify 
it with that of the harmonic ones, recalling the fact that the harmonic perturbation 
system is formally equivalent to the limit A — > of the odd system.) 

We define the stability measures by means of zero mode normalization by 

(2.35) M A = Q A /Q { °\ (A=E, O, H, or T) 

where we think of Qq = Qjj as noted above. From the asymptotic solutions for 
Qa we immediately obtain the following behaviors for our stability measures. 

Theorem 2.5. The stability measures by means of zero mode normalization for 
the vacuum orthogonal Bianchi III spacetime have the following asymptotic decaying 
or growing rates: 

Generic (m ^Q) case: 



f 0(t) (A=E), 
(2.36) \M A \ = I 0(1) (A=0 or H), 

[0(t-i) (A~T). 
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( The rates are common for all eigenvalues ( except m — 0), but depend on the kind 
of perturbations.) 

U{1)- symmetric (m = 0) case: 

f 0(t-Va) (A 2 > |) 

(2.37) |M A |= lOpft) _ (A 2 = I) 

[ 0(t -i/2+Vb^) (A 2 <i) 

for A=E or O. (The rates are common for the even and odd kinds, but there exists 
a critical eigenvalue X 2 — 1/4 for each kind.) 

Remark 2.6. One can replace (without any modifications) the time variable t 
to the proper time r in the rate formulas of the above theorem. 

Note that we have three kinds of behaviors for the generic case; The tt per- 
turbations are decaying, so we can say that the solution is stable against this kind 
of perturbations. The odd and harmonic perturbations approach a constant, so in 
this case the solution is quasi-stable. Finally, the even perturbations are growing, 
so the solution is unstable against the even perturbations. 

On the other hand all the L/(l)-symmetric perturbations are decaying, so we can 
say that the solution is stable against the £/(l)-symmetric perturbations. Although 
this means that the [/(l)-symmetric perturbations play no major role in the whole 
perturbations of our system, the existence of the critical base eigenvalue A 2 = 1/4 
is worth pointing out. A similar existence of critical value has also been reported 
in a nonlinear (but C/(l)-symmetric) analysis ]CBM . 

3. Field equations on a closed Bianchi II solution 

3.1. The background solution. The Bianchi II (or Nil) group G\i is the 
3-dimensional Lie group which consists of all real upper triangular matrices of the 
/l x z\ 

form 1 y with the usual multiplication rule for matrices. To save space 

\o ij 

let us write an element x 6 G\\ in (transposed) vector form x = (x, y, z). Then 
the multiplication rule is given by (a, b, c) o (x, y, z) = (a + x,b + y,c + z + ay) for 
another a = (a, 6, c) £ G\\. The group G\\ acts on our universal covering spatial 
manifold M = R 3 from the left, identifying the group manifold G\\ with M . The 
Bianchi II algebra is generated by the following three vectors £/ (J = 1 ~ 3) 

i\ ^9 d ^ d ^ d 

(3.1) = -K- + y-^-, ?2 = -5-) 6 = -5-- 

ox oz ay oz 

which serve as Killing vectors on our background space or spacetime. The invariant 
vectors xi an( i their dual l-forms a 1 are given by 

d d d d 

(3.2) ^ dx' ^ 2 dy dz' ^ 3 dz' 

c 1 = dx, a 2 — dy, a 3 — dz — xdy. 

Like the Bianchi type III solution the Bianchi type II vacuum solution can be 
expressed as 

(3.3) (M x R, g ab ) = (M x R, g a6 )/r, 
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using the universal cover solution (M x R, g ao ) and a covering map T acting on it. 
The spatial manifold M here is, for defmiteness, specified to be the "circle bundle 
over the torus (T 2 ) with Euler number e = 1." (See, e.g., |HS| .) The fundamental 
group can be represented in the standard notation as 

(3.4) m(M) = (.91,32,53; [51,52] = 53, [51,53] = 1, [52,53] = 1) , 

where the brackets stand for commutators, [a, b] = aba~ l b~ l . The vacuum metric 
g a b on the universal cover is given below. T C Gn is a discrete subgroup of the 
Bianchi II group Gn, which acts on the universal covering spacetime (M x R, g ao ) 
from the left in the way that keeps the natural homogeneous slicings. Some details 
for the compactification are given in the next subsection. The metric g a b onMxR 
is the induced metric from g ao . 

The exact vacuum solution g ao is given TATT] by 

(3.5) ds 2 = -N 2 (t)dt 2 + q^tXa 1 ) 2 + q 2 (t)(a 2 ) 2 + q 3 (t)(a 3 ) 2 , 
where 

(3.6) N 2 = 1 + (3 2 t Ap3 , q x = t 2pi N 2 , q 2 = t 2p2 N 2 , q 3 = 16p 2 3 P 2 t 2p3 /N 2 . 
Pi(i = 1, 2, 3) and /3 are constant parameters such that (3 > 0, p 3 ^ 0, and 

(3.7) E Pi = 2> 2 - 1. 

When pi = P2, or equivalently when qi{t) = q2(t), the solution is said to be LRS 
as in the Bianchi III case. While there seem two possible such cases (pi,p 2 ,P3) = 
(0, 0, 1) and (2/3, 2/3, —1/3), these two solutions represent the same one-parameter 
solution. When we consider a LRS solution, we may want to take (pi,P2,P3) = 
(2/3, 2/3,-1/3), since the time coordinate t in this solution approaches the proper 
time at future infinity, which is more favorable for comparisons with other models. 

As shown in |TK H . T is a four-parameter group. Our spatially closed solution 
(|3.3I) therefore comprises a six parameter solution (since the universal cover has, as 
we have seen, two independent parameters). 

Notice the fact that T is a subgroup of the Bianchi II group Gn, not necessarily 
a subgroup of a larger group like in the Bianchi III case. This in particular means 
that each invariant vector \i is we U defined not only on the universal cover M but 
also on the compactified manifold M (i.e., the induced vectors tt*xi f° r the covering 
map it : M — > M are well defined on M . We omit it* for simplicity, though). This 
is because they are invariant under the action of Gn, and so are under T C Gn. We 
will see that because of this property the invariant vectors \i pl a Y a central role in 
developing calculus concerning mode expansions. This significance of the invariant 
vectors is common for all Bianchi class A models. 

3.2. Compactifications and eigenfunctions. Consider an arbitrary locally 
homogeneous closed 3- manifold (M,q ao ) with tt\{M) given by (|3.4|) . Such a Rie- 
mannian manifold is isometric to the quotient N 5 /A, where A/" 3 = (R 3 , e 2a q a °^) 
is the universal cover with a metric that is conformal to the Bianchi II standard 
metric q^ = {a 1 ) 2 + (a 2 ) 2 + (c 3 ) 2 . e 2a is a constant scale factor. A C IsonxA/ 3 is 
a three parameter infinite group generated by three generators a^ 6 Gn C IsomA/" 3 
(i = 1 ~ 3), which are of the form (cf. KTH]) 



(3.8) ai = (u, 6, 0), a 2 = (0, 2ttv, 0), a 3 = (0, 0, 2-kuv). 
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We denote A = {ai,a2,ag}. The three parameters u, v, and S are called the 
Tcichmiiller parameters. 

Also, consider another Riemannian manifold (M, q a b), which can be represented 
as Af 3 /A, where A = {a 2 , a 3 }. (M, q a b) is a covering of (M, q a b). M is homeomor- 
phic to the direct product T 2 x R, where T 2 is the two-torus. 

Note that the (scalar) Laplacian A with respect to the standard metric 
can be expressed as A = (yi) 2 + (X2) 2 + {xz) 2 ■ This operator is apparently well 
defined on both M and M, since so is each \i- It is straightforward to confirm that 
A commutes with \3- We can therefore diagonalize eigenfunctions with respect to 
both A and X3- Consider another operator £2 = 9/dy, which we can find commutes 
with both A and X3> so one can diagonalize the eigenfunctions with respect to £ 2 , 
too. This operator however is not well defined on M, but on M, Because of this 
fact, it is convenient to first consider eigenfunctions on M, and then construct those 
on M, taking linear combinations of the eigenfunctions on M. 

Let us define eigenvalues fi and v for the operators X3 and £2 by the following 
relations 

(3.9) X3<t> = i^(f>, £,2<f> = ivcj), 

where <f> is an eigenfunction on M . Also, we define A by A(f> — —X 2 4>. From the 
boundary condition (j>(Tx) — ^(x), it is found (f> — X(x)e lfJ,z e wy , where 

u = m/(uv), m = 0, ±1, • • • , ±00, 
v = n/v, n = U, ±1, • • • , ±00. 

We call m (or ji) the fiber eigenvalue, n (or v) the auxiliary eigenvalue, X the 
total eigenvalue. The function X(x) must satisfy the following equation (that is 
equivalent to the harmonic oscillator Schrodinger equation): 

(3.11) X" + (A 2 - /i 2 - (jxx + vf)X = 0. 

This equation tells us in particular that like the Bianchi III case, one must deal 
with the generic (m 7^ 0) case and the {7(l)-symmetric (m = 0) case separately. 

Let us focus on the generic (to ^ 0) case. One can define the eigenfunctions on 
M as 



(3.12) ; , m .„(x) = Dl (±J — (—x + -))e'-V^, 

y \m\ uv v 

where -D;(C) = e~3? Hi(Q is the parabolic cylinder function defined using the 
Hermite function Hi(Q = (— l) l e^ ^re - ^ ■ The Di factor in the above equation 
satisfies Ea. H3.11|) . The plus and minus signs, respectively, correspond to to > 
and to < cases. The index I takes non-negative integers 

(3.13) Z = 0, 1, • • • ,00. 

This is related to the total eigenvalue through A 2 = (21 + + 1). 

The eigenfunctions on M are, as stated, expressed as infinite sums of these 
eigenfunctions on M. 

Theorem 3.1. Consider the generic (m ^ 0) modes. Using the eigenfunctions 
0;,m,n( x ) on the covering (M — T 2 x R, q a b), the eigenfunctions on the closed 
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manifold (M, q a b) are represented as the infinite sum 

oo 

(3.14) <&„(x)= ]T e w ("° fe+m ^^, m , no+rofe (x), 

k— — oo 

where I = 0, 1, • • • , oo, \m\ = 1, 2, • • • , oo, no = 0, 1, • • • , |m| — 1. 

Proof. Since the summand is well defined on M, which means it is invariant 
under the action of a.2 and a3, all one has to confirm is the invariance of the left 
hand side under the action of ai . It is easy to see the following transformation rule 



St, m ,„(ai x ) = A(±. (z + u) + 2) )e i&(-+«v) e iS(v+*) 



\m\ uv v' 



(3 ' 15) = e*5'A(±4fe(-« + ^S'e^ 

y |m| to v 

From this, one can confirm the invariance 

oo 

(pl,m,no l a l ° X J — 2^ W,m,no+m(fc+l) W 

/c— — oo 

(3 - 16) = E e ^ ( - (fe+1)+m ^ ) 0, m , n „ +m( , + i)(x) 

fe=— oo 

(The last equality follows from the replacement k — > fc — 1.) □ 

Note that as a result of the complete compactiheation, the index no for the 
eigenfunction on M is bounded by m. 

These eigenfunctions have the following properties. 

Theorem 3.2. For the eigenfunctions for the generic (m =/= 0) modes, the 
following is fulfilled (indices m and no are for simplicity suppressed) : 

(3.17) Am" 

X 2<t>i =±ty J | L [&4-i + tyj_i], 

where fi is defined in Eq. ]3.1(J\) and the plus and minus signs correspond, respec- 
tively, to m > and m < 0. 

PROOF. Use-D z '(C) = — |(Z?; + i(C)-?A-i(C))- The remaining task is a straight- 
forward computation. □ 

Remark 3.3. These relations can also be derived in a purely algebraic manner. 1 



4n fact, defining Ai = xi + *X2 and A2 = Xi — iX2, one can check commutation relations 
(3.18) [A,Ai] = 2i.4iX3, [A, As] = -2iA 2 X3, 
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3.3. Scalar field. As a direct application it is of great interest to see the 
scalar field equation on our Bianchi II background. For simplicity, let us consider 
a background solution where r = A, i.e., T is specified by three parameters as in 
the form of A. We use the symbols u, v and S as in the A. These parameters are 
constants throughout the spacetime. (This is not to say the Teichmuller parameters 
for a Cauchy surface are constants of motion. [TKH]) In such a case, we can use 
the eigenfunctions obtained in the previous section without any modifications. We 
expand the scalar field <E> as follows: 

(3.21) = a J,fn,no(*)0Z,m,n o (x)- 

Z,m,no 

Since the operators \i do not change m and no when acting on (f>i t m,n , we suppress 
these subscripts and write simply a; and <pi, as far as no confusion occurs. 
The scalar field equation = g ab V a Vb < I > on our background reduces to 

(3-22) = -±={^—gN-^)+ (q^ixif + feW + ^W)*, 

where a dot stands for d/dt. Using the background spacetime solution for Bianchi 
II and the y-relations, we obtain the following wave equations: 

(3.23) ai + hi + nf{t) ai = |(t- 2pi - t- 2p2 )[a,- 2 + (l + 2)(l + l)a l+2 ], 
where 

(3.24) /(f) = ^-V 2pi + r 2 ^) + 1 * (l + ph^ft-^. 

(In Ea. H3.23fl . a;<o should be considered as zero.) 

A notable property of these equations is that they comprise a set of infinitely 
coupled equations, unless the solution is LRS. These couplings between different 
modes are an unavoidable feature for the field equations on a non-LRS spacetime 
solution, because for those backgrounds the mode functions defined with respect 
to a standard metric are no longer eigenfunctions with respect to the Laplacian 
for each spatial section. That is, those eigenfunctions defined on the standard Ric- 
mannian manifold do not remain eigenfunctions in the course of the anisotropic ex- 
pansion. One could define another "time-dependent" eigenfunctions to get around 
this difficulty, but in that case time-derivatives of the eigenfunctions appear in the 
field equation, so in any case the complexities caused by anisotropic expansion are 
unavoidable. 

The above wave equation gives a closed equation when the background is LRS. 
In this case we can again find asymptotic solutions: 



which tell us, e.g., that 

(3.19) AAi<p = ([A, Ai] + Ai A)<£ = (2iA lX3 + A\A)<p = -(2fi+ A 2 )Ai0. 

This means that Ai<f> is an eigenfunction for eigenvalue A' 2 = X 2 + 2fi, and the operator A\ acts as 
an increment (decrement) operator for m > (m < 0). Similarly, A2 is a decrement (increment) 
operator with A' 2 = A 2 — 2fi for m > (m < 0). Using these relations and following a similar 
argument for, e.g., the SO(3) case to determine the ground state, one reaches at the following 
relations 

M,m = -v' 2 / 1 f !, [+l,m, M<t>l ,m = V 2 /"^!-! , m , for m > 

(3.20) ' 

Al<pi, m = V 2 M*0I-l,m, M4>l,m = — V 2 H 4>l+l,m., for m < 

which are equivalent to the relations claimed. 
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Theorem 3.4. On the LRS Bianchi II vacuum solution withpi = p2 = 2/3 and 
P3 = —1/3, the scalar field equation I3.2ty for the generic mode has the following 
fundamental solutions as t — > oo : 



Proof. As remarked in the previous section, it is a suitable choice of new 
time variable (denoted s) that is essential to have an asymptotic solution. An 
appropriate choice for the present equation is 



One can show the claim following the same procedure as in THEOREM l2.2l or l2.3l □ 

3.4. Vectors and electromagnetic field. Let us consider how we can per- 
form separation of variable for vector fields. An advantage that our Bianchi II 
model has and is common for other class A spatially closed Bianchi models is that 
in principle, we do not have to define suitable eigenvectors for this purpose. This is 
because there exists a set of well-defined independent three invariant vectors \i on 
the compactified manifold M. Any vectors on M can be expanded with respect to 
these invariant vectors (i.e., one can consider the components with respect to the 
frame {xi}) in a group invariant way and so all one has to do is to mode-expand 
the components using the already defined eigenfunctions <fii,m,n - However, this 
straightforward expansion is found not to be very useful for a successful analysis. 

A nicer way is to use an analogy with Regge- Wheeler's spherically symmetric 
case or one with Bianchi III hyperbolically symmetric case, where one can take 
advantage of the spherically or hyperbolically symmetric planes. In those cases, 
one is able to define so-called the 'even' vectors and 'odd' vectors on the symmetric 
plane (as seen in the previous section for the Bianchi III case), and the even and 
odd part equations for a linear vector field equation decouple from each other. 

Although there do not exist similar symmetric planes for the Bianchi II case, 
we can, instead, consider the "plane field" spanned by xi and X2- This plane field 
is not integrable, because the commutator [xiiXa] = X3 is independent from xi 
and X2- Nevertheless, the plane field can play a special role, because the metric 
can be LRS only when qi = q%. 

We define the area two-form e of the plane field by e a b — 2cr 1 [ a <7 2 f ) ]. We 
may define the even vectors by S a = (xi</ , ;)°' 1 a + (X20z)c 2 a , and odd vectors by 
V a = ie a b db4>i = i((x24>i)cr 1 a - {xi4>i)o~ 2 a) ■ (To raise an index for e a b we use the 
standard Bianchi II metric.) We also need orthogonal vectors — (X3<l ) i) ~ 3 'a and 
time-like vectors Sj = N~ 1 (f>i(dt) a , which may be considered to belong to the even 
part. These four kinds of vectors complete our set of eigenvectors. 

Because of the fact that the commutator of xi and X2 does not belong to the 
tangent space spanned by Xi and xi-, the even and odd parts defined this way do 
not decouple from each other, even when the background is LRS. Nonetheless, the 
use of even and odd vectors is found useful, especially when the background is LRS, 
since in that case each mode is found to decouple from the others. (Do not confuse 
couplings between the even and odd parts and ones between modes.) When the 
spacetime is not LRS, however, we will have couplings between different modes due 
to the anisotropic expansion. 



(3.25) 




(3.26) 
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Using these mode vectors, we can expand the vector potential for an electro- 
magnetic field as follows: 

(3.27) A a = J2 a (t)Sl + a 1 (t)S a + a 2 {t)V a + a 3 (t)S^. 

l,m,no 

The four kinds of functions oti{t) (i = ~ 3) serve as the field variables. Let us 
consider the LRS case below, on which background, as mentioned above, couplings 
between different modes vanish. 

The quantities we are most interested in are the (t/(l)-)gauge-invariant vari- 
ables, which can be easily found by looking at components of the field strength 
F a b = d a A — d A a . We obtain the following four independent gauge invariant 
variables: 

(3.28) Qi = a>i — a 3 , Q 2 = a 2 , P\ = di — a a , P 2 = d 2 . 

Although function P% = a 3 — ao is also gauge invariant, it is found that it can 
be (consistently) solved with the others, due to the constraint part of Maxwell's 
equation. The evolution equations for the LRS background are found as follows: 

1 93, 



(3.29) 




l)Pl+P 2 ), Q2=P2, 

* N2 n 

- \i — Qi, 

93 

N 2 N 2 

- M — Q2 - M — (Qi + (21 + 1)Q 2 ). 

93 9i 



To tackle the gravitational perturbation equation we need to generalize our 
framework so that it is applicable to systems with symmetric tensors. Again, thanks 
to the well-defined invariant vectors xii separations of variable for the tensors are 
straightforward in principle. However, because of the complexity of the linearized 
Einstein equation, to make our analysis successful we will need to find a well- 
arranged set of base mode tensors like the even and odd vectors used above. 



References 

[CBM] Y. Choquct-Bruhat, and V. Moncrief, Future global in time Einsteinian spacetimes with 

U(l) isometry group Ann. Henri Poincare 2 (2001) 1007—1064; Non linear stability of 

einsteinian spacetimes with U(l) isometry group, preprint ( gr-qc/0302021 1. 
[CH] L. Cesari, Asymptotic behavior and stability problems in ordinary differential equations, 

Springer, 3rd ed. (1971); P. Hartman, Ordinary differential equations, John Wiley (1964). 
[EM] G.F.R. Ellis, and M.A.H. MacCallum, A class of homogeneous cosmological models, Com- 

mun. Math. Phys. 12 (1969) 108-141. 
[FM] A.E. Fischer, and V. Moncrief, The reduced Einstein equations and the conformal volume 

collapse of 3-manifolds, Class. Quantum Grav. 18 (2001) 4493—4515; and reference therein. 
[HS] J. Hempel, 3-manifolds, Ann. of Math. Studies 86, Princeton University Press (1976); 

P. Scott, The geometry of 3-manifolds, Bull. London Math. Soc. 15 (1983) 401-487. 
[KTH] T. Koike, M. Tanimoto, and A. Hosoya, Compact Homogeneous Universes, J. Math. Phys. 

35 (1994) 4855-4888. 

[TAU] A. H. Taub, Empty space-times admitting a three-parameter group of motions, Ann. Math. 
53 (1952) 472-490. 

[TKH] M. Tanimoto, T. Koike, and A. Hosoya, Dynamics of compact homogeneous universes, J. 
Math. Phys. 38 (1997) 350-368; Hamiltonian structures for compact homogeneous uni- 
verses, J. Math. Phys. 38 (1997) 6560-6577. 



16 



MASAYUKI TANIMOTO 



[TMY] M. Tanimoto, V. Moncrief, and K. Yasuno, Perturbations of spatially closed Bianchi III 

spacetimes, Class. Quantum Grav. 20 (2003) 1879-1927. 
[WE] J. Wainwright, and G.F.R. Ellis, Dynamical Systems in Cosmology, Cambridge University 

Press (1997). 

[YT] K. Yasuno, and M. Tanimoto, Asymptotic behavior of linearly perturbed spatially closed 

Bianchi III vacuum solutions, in preparation. 
[YO] J.W. York Jr., Conformally invariant orthogonal decomposition of symmetric tensors on 

Riemannian manifolds and the initial-value problem of general relativity, J. Math. Phys. 

14 (1973) 456-464. 



Errata to the published version 

This article will be published in Contemporary Mathematics of American Mathe- 
matical Society. The following is corrections to that published version. (Unfortunately, 
authors to this publication are not given a chance to make corrections after the final 
submission.) All the corrections have already been incorporated in the present electronic 
version, so there is no need to make any change for this article. 

1. The line below Eq.(3.14) should read 

"where I = 0, 1, • • • , oo, |m| = 1, 2, ■ ■ ■ , oo, n = 0, 1, • • • , |m| — 1." 

2. Eq.(2.2) and Eq.(3.1) should be swapped. 
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